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We study the superfluid state of two-species heteronuclear Fermi gases with isotropic contact and anisotropic 
long-range dipolar interactions. By explicitly taking account of Fock exchange contribution, we derive self- 
consistent equations describing the pairing states in the system. Exploiting the symmetry of the system, we 
developed an efficient way of solving the self-consistent equations by exploiting the symmetries. We find that 
the temperature tends to increase the anisotropy of the pairing state, which is rather counterintuitive. We study 
the anisotropic properties of the system by examining the angular dependence of the number density distribution, 
the excitation spectrum and the pair correlation function. The competing effects of the contact interaction and 
the dipolar interaction upon the anisotropy are revealed. We derive and compute the superfluid mass density 
Pij for the system. Astonishingly, we find that p zz becomes negative above some certain temperature T*(T < 
T c ), signaling some instability of the system. This suggests that the elusive FFLO state may be observed in 
experiments, due to an anisotropic state with a spontaneously generated superflow. 

PACS numbers: 03.75.Ss, 05.30.Fk, 67.85.Lm, 74.20.Fg 



I. INTRODUCTION 

The recent experimental realization and coherent control of 
high phase-space density quantum gas of the polar molecules 
40 K 87 Rb HHj] has stimulated great interests in studying 
the effect of anisotropic long-range dipole-dipole interactions. 
Furthermore, the large dipole moment of the polar molecules 
can be tuned using an external electric field; This has a range 
of applications such as the control of ultracold chemical re- 
actions y[], the design of a platform for quantum information 
processing |4|, dD and the realization of novel quantum many- 
body systems laLZI]- 

For the fermionic polar molecules, dipolar interactions has 
led to various new phenomena, including fractional quantum 
Hall effect [8], Wigner crystalization [9], nematic phases 1 1011. 
interlayer superfluidity lUlll and density wave patterns IU2I1 . 
Two fundamental properties of dipolar Fermi gases are super - 
fluid pairing fl3l4l8ll and Fermi surface deformation IU9T - |22|] , 
originating from the partially attractive nature of the dipolar 
interaction and anisotropic Fock exchange interaction. For 
dipolar Fermi gases with two hyperfine states, one can tune 
not only the dipole-dipole interaction by a fast rotating ori- 
enting field [23], but also the s-wave interspecies interaction 
via a Feshbach resonance. Therefore, one naturally expects 
rich physics will emerge as a result of the interplay of the 
anisotropic long-range dipole interaction and short-range s- 
wave interaction. 

There have been intense theoretical efforts on the proper- 
ties of the normal phase at zero temperature, including zero 
sound propagation 02411 , Fermi liquid description I125ll . and ex- 
pansion dynamics 11201 l26j - l28ll . At finite temperature, there 
have been several theoretical efforts. Endo Il29ll et al. devel- 
oped a variational approach, which allowed them to assess the 
effect of temperature on the position- and momentum-space 
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distortion. Zhang and Yi [30] explored the system deforma- 
tion and stability by using a Monte Carlo method. Baillie and 
Blakie [22] studied thermodyamics and correlation properties 
of the system by employing a semiclassical self-consistent 
Hartree-Fock calculation. 

In this paper, we consider BCS pairing of two-species dipo- 
lar Fermi gases at finite temperature, aiming to investigate 
thermodynamics and anisotropic properties in the superfluid 
state. To fully assess the effect of Fock exchange term, it is 
essential to explicitly take account of the Fock exchange term 
in a self-consistent way. The content of the paper is organized 
as follows. In Sec. II, we present our model Hamiltonian for 
a homogeneous two-species dipolar Fermi gas and develop a 
mean-field description for the superfluid state at finite temper- 
ature. In Sec. Ill, We present theoretical studies and numerical 
results on the anisotropic properties and the thermodynamics 
of the system . In Sec. IV, we derive and compute superfluid 
mass density, then explore the consequences. Finally, we con- 
clude in Sec. V. 



II. MODEL HAMILTONIAN AND MEAN-FffiLD 
DESCRIPTION 



We consider a homogeneous gas of two species of 
fermionic heteronuclear molecules a =t and |. For sim- 
plicity, we further assume that each species has the same 
mass, density and dipole moment. The electric dipoles of the 
molecules with moment d are oriented along the z axis by 
a sufficient strong external electric field such that the spin- 
independent part of the electronic dipole-dipole interaction 
becomes Vdd{q) = (4?r/3)(3 cos 2 # q — 1), with # q being the 
angle between momentum transfer q and the direction of z 
axis in which the dipoles are aligned. In addition, we assume 
that the molecules also interact via a contact interaction with 
strength g. This system is described by the following Hamil- 
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where /i is the chemical potential, n is the total number den- 
sity, V is the volume, and = k 2 /2m (where we have set 
h = 1). The interaction potential V aa ' (q) = gS a .- a i +Vdd(q) 
contains both dipole-dipole and contact interactions. Antici- 
pating the importance of the Fock exchange term, we decou- 
ple the interaction in all three channels [31]: direct channel, 
exchange channel and Cooper channel, resulting in the fol- 
lowing effective mean-field Hamiltonian 
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Here £ k o- = £k — n + gn/2 + Sk CT , with self-consistent mean- 
fields defined as 
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Some comments are in order: The contact interaction affects 
the single-particle spectrum by shifting the chemical potential, 
which may be redefined as fi = // — gn/2. The self-energy 
E k(T encodes the anisotropic dipolar contribution from Fock 
exchange term to the dressed single-particle spectrum, which 
justifies our treatment. In addition, both parts of the interac- 
tion contributes to the pairing field. 

The Hamiltonian (f2j) is diagonalized by invoking the Bo- 
goliubov transformation 03211 . We obtain self-consistent equa- 
tion for the self-energy, the order parameter and the number 
density 
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where E ka = v^L + J2<r> |A CT(T '(k)| 2 = E k is the quasi- 
particle spectrum and /3 = 1/fcfjT is the inverse tempera- 
ture. Equation © formally diverges. For the contact inter- 
action, one can eliminate the interaction strength g in favor 
of the s-wave scattering length a s using 1/g — m/(Aira s ) — 
1/V^ k (2e k ) _1 . The dipolar interaction can be regularized 
by replacing the bare interaction V aa ' (k — p) with the vertex 



function r CTCr '(k — p) as explained in Refs. lfl3i[33ll . To first 
order in Born approximation, the gap equations become 
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The above equation, together with Eqs. (O and (Q, comprises 
a complete description of the dipolar Fermi gas and needs to 
be solve self-consistently. 

Due to the symmetry of the interaction potential, the mo- 
mentum distribution, order parameter and self energy possess 
azimuthal symmetry in thermal equilibrium, as could be seen 
from the self-consistent equations by integrating out the az- 
imuthal degrees of freedom k . Therefore physical quantities 
such as n k , A (T(T /(k) and S kcr are only functions of (fc,#k)- 
Numerically, we parameterized these functions by two dimen- 
sional grids living on the domain f2 = [0, k c ] x [0, ir] (where 
k c is the momentum cutoff). In this study, we focus on the 
spin singlet BCS pairing, and suppress the spin index for or- 
der parameter by writing A^(k) = A(k). The singlet or- 
der parameter possesses inversion symmetry as well as az- 
imuthal symmetry A(/c, k ) = A(fc,7r — # k ). It is easily to 
verify that n k and E k share the same symmetry. Thus the 
results presented below is limited to < #k < n/2. In 
our calculation, we parameterize the dipolar interaction by 
the dimensionless coupling parameter Cdd = md 2 {n a ) 1 / 3 ' ', 
where n a = k F /(Qir 2 ) = n/2, and the Fermi energy is 
Ep = k 2 F /2m. To make our mean-field study qualitatively 
valid and quantitatively reliable, we will restrict ourself to 
studying weak-coupling BCS regime. 



III. ANISOTROPIC SINGLET SUPERFLUID AT FINITE 
TEMPERATURE 

In the normal phase at zero temperature, the Fermi sur- 
face of the dipolar gas has an ellipsoid shape lfl9ll : n k = 
Q(kp — a 2 k 2 — k 2 /a — k 2 /a), where a is the deforma- 
tion parameter. In a previous study 111 811 . we generalized 
the definition of the deformation parameter a to measure the 
anisotropy of the single particle momentum distribution with 
a similar strategy. The angular distribution ng(9-^) can be 
obtained by integrating out the magnitude and azimuthal an- 
gle of the momentum n#((9k) = J n(k)k 2 dk/(2ir 2 ). Since 
ng(Q)/ng{-K /2) = k\jk\ and a 2 k 2 = k 2 /a,we deduce that 
a = [rig (vr/2)/n e (0)] 2 / 9 . In Fig. 5] the deformation param- 
eter a is plotted as a function of the temperature. For a fixed 
for s-wave coupling strength 1/kpa — —0.6, we show, in the 
panel (a), the deformation parameter a for three typical dipo- 
lar coupling strength Cdd — 0.6, 0.8 and 1. In each curve, 
the deformation parameter a develops a minimum at a critical 
temperature T c separating the superfluid state from the nor- 
mal state. The calculated critical temperature T c is 0.175Tp, 
0.166T F and 0.155T F for C dd = 0.6, 0.8 and 1, respec- 
tively. In the normal phase where (T > T c ), a increases as 
the temperature increases. This is expected as increasing the 
temperature tends to make the system more isotropic by re- 
ducing the effects of anisotropic interaction. In the superfluid 
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FIG. 1 . (Color online) The temperature dependence of the deforma- 
tion parameter a. (a) At s-wave coupling strength 1/fcf-a = —0.6 
and for different dipolar coupling strength Cdd', (b) At dipolar cou- 
pling strength Cdd = 0.6 and for different s-wave coupling strength 
1/kpa. The minimum point in each curve corresponds to the posi- 
tion of phase transition between the superfluid state and the normal 
state. Tf is the Fermi temperature. 



state, the deformation parameter a decreases as the tempera- 
ture increases, which is rather counterintuitive. This may be 
understood in the following way: increasing the temperature 
tends to weaken the BCS singlet pairing state which favors 
an isotropic momentum distribution. In both the normal state 
and the superfluid state, for fixed s-wave coupling, the larger 
the dipolar coupling strength Cdd, the smaller the deforma- 
tion parameter a. For a fixed Cdd = 0.6, as shown in the 
panel (b), the deformation parameter decreases as the temper- 
ature increases for all three typical s-wave interaction strength 
1/kpa — —0.6, —0.7 and —0.8. The critical temperature 
T c as well as the deformation parameter increases as one in- 
creases the s-wave coupling strength. In the normal phase, 
increasing the s-wave interaction strength has no effect on the 
deformation parameter, as the pairing channel closes and the 
effective chemical potential remains the same so as to fixed 
the average particle number. 

To further reveal the anisotropic nature of the pairing states, 
we investigate the pair correlation function and the quasipar- 
ticle excitation spectrum. The pair correlation function can be 
evaluated as [ 34] 

Ctj.(ki,k 2 ) = < nu^nu^i > - < rikit >< ™k 2 i > 



tanh 




As shown in Fig. [2^, the quasiparticle excitation spectrum 



FIG. 2. (Color online) (a) The quasiparticle excitation spectrum 
and (b) the pair correlation function Ck along various polar angles as 
a function of momentum magnitude for the s-wave coupling strength 
l/k F a — —0.6 and the dipolar coupling strength Cdd = 1 at zero 
temperature. 



Ek starts at the same value from zero momentum and evolves 
anisotropically in momentum space. Along any polar angle, 
i?k decreases at small momentum magnitude and increases 
at high momentum magnitude; it develops a minimum which 
is greater than zero, signaling that it is a gapped superfluid. 
The locus of the minimum varies from a small momentum 
magnitude to a large momentum magnitude as the polar angle 
6k decreases from ir/2 to 0, indicating a distortion along the 
z axis. The pair correlation function Ck is shown in Fig. \2)p. 
It is interesting to notice that for a fixed polar angle, there 
exists some typical momentum magnitude at which the pair 
correlation function Ck achieves a maximum. Remarkably, at 
this typical momentum, the quasiparticle excitation spectrum 
i?k develops a minimum. As the polar angle changes from 
to 7r/2, the momentum magnitude at which the maximum 
pairing occurs decreases. At a large momentum magnitude, 
the pair correlation decays rapidly with 1 /fc 4 asymptotically. 

To assess the temperature effect on the anisotropic prop- 
erties of the system, we study momentum space distribu- 
tion of the number density. The plot in panel (a) and (b) 
of Fig. [3] corresponds to radial momentum distribution of the 
number density along the polar (6*k = 0) and axial direction 
(#k = tt/2), respectively. At 6>k = 0, the number occupation 
at any given momentum magnitude increases slightly as the 
temperature is raised from T/T F = 0.05 to T/T F = 0.15. 
In contrast to the polar direction, the situation at #k = tt/2 
looks quite different: for k < kp, the number occupancy 
for T/Tp = 0.05 and T/T F = 0.15 coincides at two typ- 
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FIG. 3. (Color online) The radial number occupancy ne(#k) at 
s-wave coupling strength 1/fcf-a = —0.6 and dipolar coupling 
strength Cdd ~ 1.0 for (a) the polar direction 8^ = and (b) the ax- 
ial direction 9^ — n/2. The solid blue curves are for T/Tf = 0.05, 
and the dash red curves are for T/Tf = 0.15. Tf = Ef /ks is the 
Fermi temperature. 



FIG. 4. (Color online) The entropy S/(NkB) as a function of 
the temperature, (a) At s-wave coupling strength = —0.6 

and for different dipolar coupling strength Cdd', (b) At dipolar cou- 
pling strength Cdd = 0.6 and for different s-wave coupling strength 
1/feFfl. The linear (high temperature) part of the curves corresponds 
to the normal phase, characteristic of Fermi liquid behavior. 



ical momentum magnitudes; For momentum magnitude lies 
in the range of0.5fci? < k < 0.6k F , the number occupancy 
for lower temperature exceeds that of higher temperature; the 
trend appears to be reversed for 0.6k F < k < 1.0k F . Judging 
from the cases of 8^ — and 8^ = 0, one can largely de- 
duce that as the temperature increases, the number occupancy 
transfers from the axial direction to the polar direction, ren- 
dering the system more anisotropic. This is consistent with 
the temperature dependence of the deformation parameter a 
shown in Fig.Q] 

The Entropy provides useful insights into the thermody- 
namics of the superfluid pairing states. At mean-field level, 
the system may be regarded as a collection of noninteracting 
quasiparticles with dispersion E ka , obeying Fermi statistics. 
Thus, the entropy of the system is evaluated as 

S = -k B ^ ["^(^kcr) lnn F (S ko .)] 

-k B K 1 - M^ka)) In (1 - n F {E ka ))} 

k<r 

= k B [/3E ka n F {E k „) + ln(l + e^-)] , (9) 

where nf is the Fermi function defined by n F (E ka ) — 
1/[1 + exp (/3-Eko-)]- In Fig- the entropy is plotted as a 
function of the temperature for fixed s-wave coupling strength 
l/k F a = —0.6 and for different dipolar coupling strength. 



As seen in the figure, for every dipolar coupling strength Cdd 
there exists a critical temperature T c separating the normal 
phase from the ordered pairing phase. The entropy in the nor- 
mal phase has characteristic linear temperature dependence, 
a signature of the Fermi liquid behavior. We find that, in the 
superfluid phase, a larger dipolar interaction strength leads to 
a higher entropy, which is opposite to what is found in a po- 
larized normal phase ll22l l30ll . As a measure of the degree of 
disorder, a higher entropy characterizes a less ordered state. 
This finding is consistent with the established argument lfl8ll 
that in the singlet pairing phase the anisotropic dipolar inter- 
action is a competing effect as opposed to the s-wave contact 
interaction. In Fig. HJ?, the entropy decreases as the s-wave 
coupling strength 1/kpa increases. When the system is in the 
normal phase, increasing the s-wave interaction has no effect 
on the entropy. This is due to the fact that, in the normal state 
the pairing channel is closed, and the effective chemical po- 
tential Ji = fj, — gn/2 is the same so as to fix the average 
particle number. 



IV. SUPERFLUID DENSITY 

The superfluid density is a fundamental quantity describing 
the response to a rotation as well as in two-fluid collisional hy- 
drodynamics. In the following, we will obtain the superfluid 
density using the thermodynamic potential in the presence of 
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FIG. 5. (Color online) The superfluid density tensor pij as a function 
of temperature for s-wave coupling strength 1/kpa = —0.6 and 
for various dipolar coupling strength Cdd' (a) pxx, (b) p zz and (c) 

p = {pxx + Pyy + Pzz) /3. 



a superfluid flow. To impose a current, one applies a "phase 
twist" to the order parameter A(x): A(x) — > A(x)e lC *' r . 
The superfluid velocity v s associated with this imposed phase 
twist is v s = Q/M, and M is the Cooper pair mass M = 2m. 
The grand potential with superfluid velocity v s could be writ- 
ten as B5[ 13611 



fi(v.) = fi(0) +j s • v s + -ftj(v«)i(v.)j + •• •, (10) 
0(v s ) - n(0) = -^(TrlnG- 1 - Tr InG -1 ). (11) 

In the above, j s is the superfluid current, pij is the super- 
fluid (mass) density tensor, the trace operator Tr runs over 
both spin space and space-time coordinate x — (r, r), and 
G -1 (x, x') is the inverse of the Green's function in the Nambu 
space 



G~\x,x') = 5{x-x') x 

f-3r-£+M-£(z) 



A*(x) ,, T — 
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To bring Tr In G s 1 into a more convenient form, we intro- 



duce a unitary matrix 
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Employing Tr In G^ 1 = Trln [WG^U] = Tr In CT 1 , we 
obtain 



G _1 (a,a/) = 5(x-a;') x 
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A(z) 

_Q\2 



(i + S(x)/ 
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Transforming Eqs. (TlZt and (fT4l to the momentum-frequency 
representation, we obtain 

G- 1 (k,iw n ) = G- 1 -±-S.®i-QL®a z . (15) 
2m 8m 

Therefore, one can perform the derivative expansion, 



Tr InG - TrlnG" 1 = k • v s TrG 

\2 
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Combining Eqs. ([Tol l. (fTTT i and ( fTol l. we readily obtain the ex- 
pression for superfluid density tensor 
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To proceed further, one needs the Green's function G(k, iw n ) 
obtained from Eq. (fTZt 



G(k, iw T: 



E. 
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•(18) 



Substituting Eq. ( fl"8l) into Eq. ( fTTI i and carrying out the Mat- 
subara frequencies summation, we finally reach 



Pij 



mn - 



(19) 



The cylindrical symmetry implies that p xx — p yy . The aver- 
age superfluid density is defined p = (p xx + p yy + p Z2 )/3 = 
(2p xx + p 2Z )/3. 

The temperature dependence of superfluid density is shown 
in Fig.|5]for fixed s-wave interaction strength 1/kpa = —0.6 
and for different dipolar interaction strength Cdd = 0.6, 0.8 
and 1. For a given Cdd, Pxx, Pyy and p decreases as the tem- 
perature increases until the temperature reaches the transition 
temperature T c . Higher Cdd weakens the pairing states, lead- 
ing to lower p xx , p zz and p. At zero temperature, all the par- 
ticles contribute to the average superfluid density p. As seen 
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in panel (c), as the temperature goes up, p drops down, and 
reaches zero at the normal phase. This behavior is similar to 
what is expected from the picture of the two-fluid model. By 
close inspection of panel (b), one finds a peculiar behavior: for 
a given Cdd, the Pzz is negative when the temperature is higher 
than some typical temperature T* (T* < Tc), which signals 
some instability of the system. The FFLO pairing state with a 
nonzero center of mass momentum along the polar direction is 
energetically favored due to the negative p zz . The search for 
the FFLO state has lasted for more than four decades in many 
branches of physics. The observation of the FFLO phase in 
solids proved very difficult and was only achieved recently 
in the heavy fermion superconductor CeCoIn5 Q7I1 . In cold 
atom physics, very recently, evidence of the FFLO superfluid 
has been found [38] in the experiment on a one-dimensional 
Fermi gas. 

V. CONCLUSION 

To sum up, we have presented a detailed theoretical and nu- 
merical study on the singlet superfluid pairing in a two-species 
dipolar Fermi system at finite temperature. Fascinating conse- 
quences of long-range anisotropic dipolar interaction are ex- 
plored. The competing effects of isotropic s-wave contact in- 
teraction and anisotropic dipolar interaction are revealed. We 
also study finite temperature effect on the system through var- 
ious physical quantities. Importantly, We find that increasing 



the temperature tends to make the system more anisotropic. In 
particular, we derive and compute the superfluid density ten- 
sor for the system. We argue that the elusive FFLO state may 
be realized in polar molecules. The phase space where FFLO 
state is stable and the properties of the FFLO state deserves 
further investigation. 

The anisotropic properties of momentum distribution and 
excitation spectrum bear consequences for experimental ob- 
servations. The single-particle momentum distribution is 
routinely observed by time-of-flight measurements ll39l . 
Anisotropic features of the quasiparticle spectrum could be 
probed by radio frequency spectroscopy l4lll . The recent 
experimental work 113 811 paves the way to direct observation 
and characterization of FFLO pairing. With the experimen- 
tal progress toward achieving degenerate polar molecules, one 
may expect that the predictions could be verified experimen- 
tally in near future. 
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